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Abstract 

We present an exactly solvable model of a Gaussian (flexible) polymer chain 
in a quenched random medium. This is the case when the random medium 
obeys very long range quadratic correlations. The model is solved in d spatial 
dimensions using the replica method, and practically all the physical prop- 
erties of the chain can be found. In particular the difference between the 
behavior of a chain that is free to move and a chain with one end fixed is 
elucidated. The interesting finding is that a chain that is free to move in 
a quadratically correlated random potential behaves like a free chain with 
R 2 ~ L, where R is the end to end distance and L is the length of the chain, 
whereas for a chain anchored at one end R 2 ~ L 4 . The exact results are 
found to agree with an alternative numerical solution in d = 1 dimensions. 
The crossover from long ranged to short ranged correlations of the disorder is 
also explored. 

PACS number(s): 05.40-a, 75.10.Nr, 36.20.Ey, 64.60.Cn 
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I. INTRODUCTION 



The behavior of polymer chains in random media is a well studied problem [1-6] that 
has applications in diverse fields. Besides the polymers themselves this problem is directly 
related to the statistical mechanics of a quantum particle in a random potential 0, the 
behavior of flux lines in superconductors in the presence of columnar defects @,@], and the 
problem of diffusion in a random catalytic environment [Q. Despite the volume of work 
that has been done on these problems there are still many unanswered questions. Most 
of the previous work (with the exception of directed polymers [|K],|Tl|]) concentrated on 



disorder with short ranged correlations. In this paper we consider a model with long ranged 
(quadratic) correlations of the random potential that can serve as a laboratory (toy model) 
since it can be solved exactly using the replica method [T^| . Since some people are somewhat 
wary of the n —>■ limit used in replica calculations, we also solve the model numerically 
in one dimension and obtain an excellent agreement with the analytical solution. More 
importantly, the numerical solution enables us to explore the crossover from long ranged to 
short ranged correlations of the disorder and obtain a coherent picture of the behavior of a 
Gaussian chain in a random medium. 

The simplest model of a polymer chain in random media is a Gaussian (flexible) chain 
|13| in a medium of fixed random obstacles |J. In this paper we do not include a self-avoiding 
interaction. This model can be described by the Hamiltonian 



cL 

H= I du 



M ( d~R(u)\ 2 ji 2 



+ ^R 2 (u) + V(R(u)) 



(1.1) 



2 V du J 2 

where R(w) is the d dimensional position vector of a point on the polymer at arc-length 
u (0 < u < L), and where L is the contour length of the chain. The medium of random 
obstacles is described by a random potential ^(R) that is taken from a Gaussian distribution 
that satisfies 

(V(R)> = 0, (V(R)V(R0> = / ((R - R') 2 ) • (1-2) 

The harmonic term in the Hamiltonian is included to mimic the effects of finite volume. 
This is important to ensure that the model is well defined, since it turns out that certain 



equilibrium properties of the polymer diverge in the infinite volume limit (/x — > 0). The 
function / characterizes the correlations of the random potential, and will depend on the 
particular problem at hand. The parameter M is inversely proportional to (3b 2 , where 
P = (fcsT) -1 , and where b is the Khun bond step. 

Once we have defined the Hamiltonian for any chain configuration R(w)> we can write 
the partition sum (Green's function) for the set of paths of length L that go from R to R' 
as 

'R(0)=R 

All the statistical properties of the polymer will depend on the partition sum. For instance, 
we can calculate the averaged mean squared displacement of the far end of a polymer with 
one end that is fixed at the origin. This is a measure of the wandering of a tethered polymer 
immersed in a random medium. This quantity can be written as 



/-R(L)=R' 

Z(R, R'; L) — [dR(u)]exp(-/3H). (1.3) 

JR(0)=R 



/ dRR 2 Z(0, R; L) 



where the overbar stands for the average of the ratio over the realizations of the random 
potential. This average is referred to as a quenched average, as opposed to an annealed 
average, where the numerator and denominator are averaged independently. For a polymer 
with one end fixed a typical conformation in a random medium is that of a tadpole. The 
head of the polymer wanders far from the origin to find a region of favorable potential and 
then the remaining chain settles itself in that region. This is at least what is believed to 
happen when the disorder has short ranged correlations HQ . On the other hand if the chain 
is not anchored but both ends are free to move, the head to tail mean squared displacement 
is given by 



/ dRdR'(R - R) 2 Z(R, R' ; L) 



In this case the chain can move as a whole to find a favorable environment in the random 
medium. 
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In order to compute the quenched average over the random potential we apply the replica 
method. We first introduce n-copies of the system and average over the random potential 
to get 

/-Re(L)=H4 « 

Z n {{R a }, {RJ; L) = Z(R ls R' i; L) ■ ■ ■ Z(R n , R' n , L) = / X\{dR a ] exp(-(3H n ), 

(1.6) 



where 



1 f 

H n = - duY] 

-§ / d«/ ^'V/((R a H-R fe K)) 2 ). (1-7) 

The averaged equilibrium properties of the polymer can now be written in terms of the 
replicated partition sum Z n ({R a }, {R Q }; L). For instance, the mean squared displacement 



defined in Eq. (O) can be written in as 



1WTTY\- I™ J^i---^R^({Q},{R a };L) 



and similarly 



( Fl j)_ -o /n^n^ / ^n({R a },{R a }^) ' 1 j 

Thus, the averaged equilibrium properties of the polymer can be extracted from an n-body 
problem by taking the n — > limit at the end. This limit has to be taken with care, by 
solving the problem analytically for general n, before taking the limit of n —>■ 0. 

We now proceed to introduce our toy model that can be exactly solved using the replica 
method and which also lends itself to an accurate numerical solution. This is the case when 
a Gaussian polymer chain is immersed in a random medium that has very long range spatial 
correlations. In particular, we take the correlation function to be of the form 

(V(R)V(R')) = f ((R - R) 2 ) = g (1 - (R - R)7£ 2 ) , (1.10) 

where £ is chosen to be larger than the sample size, so that the correlation function is well 
defined (non-negative) over the entire sample. Since this model can be solved exactly using 



the replica method, we can compute all the important physical properties of the polymer 
chain, and then compare the exact analytical results with an alternative numerical solution 
(at d — 1). Also, this model of long range correlations is interesting in its own right in 
that it may serve as a good approximation to any correlation function / that is smooth and 
slowly decaying. Most cases investigated so far in the literature are concerned with disorder 
with short ranged correlations. 

There are many properties of the polymer chain that can be exactly computed. In 



addition to (Rj-(L)) and (Hp(L)) we will compute two other quantities. First, for a polymer 
loop of arc-length L, we will compute the quantity 



in the limit L ^> I. This is a measure of the average fluctuations of a chain segment of 
arc-length /. Since in this case the chain is not anchored, this quantity is in some respect 



similar to (R F (L)). Yet another quantity of interest is 



/ dKR 2 Z(K, R; L) 



(RqW) = y j (fltZ(R, R; £) J ' (L12) 
which has a more direct application to the related problem of a quantum particle in a random 
potential [|7| . The reason for this is that the partition sum of a polymer chain can be mapped 
to the density matrix of a quantum particle. The mapping J7|,[I|J is given by 

P -> 1/h, L -> (3h. (1.13) 

Then p(R, R'; (3) = Z(R, R'; L = /3h, (3 = 1/h) is the density matrix of a quantum particle at 
inverse temperature (3. Note that the variable u is now interpreted as the Trotter (imaginary) 



time, and M as the mass of the quantum particle. Under this mapping (Rq(L)) can be 
interpreted as the average mean squared displacement of a quantum particle in a random 
plus harmonic potential. 

The paper is organized as follows: In Sec. || we outline the exact analytical solution for 
various quantities relevant to a polymer chain. In the next section we present the details of 



the numerical approach to the problem. In Sec. [V we compare the analytical and numerical 



results and comment on the physical implications of our results. Concluding remarks are 
offered in Sec. |V|. 



II. THE ANALYTICAL SOLUTION 



We start with the case when one end point is fixed. The analytical calculation is based 
on an exact evaluation of the replicated partition sum (|1.6| ). For the correlation function / 
that we are considering the replicated Hamiltonian is 



(2.1) 



where a = g/2£ 2 , and where we have dropped the constant part of the function / since it 
only contributes an unimportant normalization factor. After expanding the quadratic term 
and simplifying the double integral we get the replicated Hamiltonian 

lT 2 "1 / l \ 2 



Now, using the Gaussian transformation 



(2.2) 



,Q 2 /2 



(27r) d / 2 



rfAe (~A 2 /2- Q .A) 



and letting 



we can write the replicated partition sum as 



Z n ({R a }, {R'J; L) 



(27T 



d/2 



dXe II / [dR a ]e-^ X \ 

a=1 Jlla{0) = Ka 



where 



H a {\) = I du 
'o 



f L , M (dR a (u)\ 2 u' . , „ 



u 



(2.3) 



(2.4) 



(2.5) 



(2.6) 
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and where // = // + 4nf3aL. The path integrals can now be evaluated directly using well 
known results for quadratic Hamiltonians. The details of the calculation are given in the 
Appendix. Once the partition sum is known we can directly evaluate the right hand side of 
Eq. (|1.8|) by taking n —>■ at the very end. The result is 



d 1 I \ fi \ Aad I 1 



(R2 (L)\ = / tanh \\ i-L + 1 , = , , (2.7) 

v TV ;/ (3 V M\x \\ M J /i 2 I cosh(^L) 1 



We can also compute the averaged mean displacement square (Rt(L)) 2 (see Appendix). 



We find that 



This implies that the displacement from the average is 



(HMD) " (Rt(D> 2 = Iv^tanh yZ L ) , (2.9) 



which is independent of disorder. We will discuss the physical implications of these results 
in a later section. 

Considering now a chain that is free to move we calculate (see Appendix for details) the 



quantity (R|(L)) using Eq. flO]). The result is 



2d sinh( v /^L) 



(Rl(L)) = 7= 7 ; V ^- M N ' v , (2.10) 

V Fl ;/ /VM^cosh^L) +1)' 1 ; 



which is independent of disorder and in the limit of /x — > behaves like dL//3M, i.e. like a 
free chain. 



The quantity (Rq(L)) can also be computed exactly from the expression 

J dUi ■ ■ ■ <fR n R 2 Z n ({R a }, {R a }; L) 



{R2q{L)) ~ -o fdR 1 ...dR n Z n ({R a h{R a }-,L) ■ (2 ' n) 
Details are given in the Appendix. However, in this case it is also possible to carry out the 
computation in an alternative way by taking advantage of the periodic boundary conditions 
of the closed loop. This provides for a further check on the result and is also included for 
instructional purposes. Using the Fourier space variables 



R a (tu) = [ duR a {u)e 

Jo 



we can write the propagator associated with (5H n as 
(3G ab (uj) = ^(R a (uj) -R b (-uj)} 



{ (Moo 2 + li + 4n(3aL) I - A/3<tL5 U)0 }J , (2. 12) 



where u is restricted to the discrete values 



2-7T 

oj m = —m, m = 0, ±1, ±2, ... . (2.13) 



After inverting the n x n matrix and taking the n — > 0, we find 



Then, using the relation 



/3G ab (a, = 0) = 5 -± + ^, (2.14) 
H /i 2 



we find that 



which implies that the only effect of the disorder is to shift the zero disorder result by a 



constant factor. Next, we compute the quantity (Rq(L)) 2 . We find that 
1 n 



a^b 



Ln(n - 1) *-f fi 



which again implies that the deviation from the average (Rq(L)) — (Rq(L)) is independent 
of disorder. 

Finally, we compute the quantity C(/), which was defined in Eq. (|1.11| ). We use, 

n 



nd 

0=1 



which for large L yields the expression 

C(l) = — != (l - exp(-/v^7M) . (2.20) 
So we find that C(l) is independent of the disorder and is the same as that of a free chain. 



III. NUMERICAL PROCEDURE 

In order to check the validity of the analytical solution we will have to numerically 
compute the quenched average of certain physical properties of the polymer. This is a 
rather computationally intensive task because of the difficulty of evaluating the partition 
sum, and also because all quantities will then have to be averaged over many realizations of 
the random potential. In this paper we will only concentrate on the case d = 1. Although 
this does not correspond to a physical polymer (d = 3) we will still be able to check the 
validity of our analytical results for the special case d = 1. In the context of the quantum 
particle in a random potential this case corresponds to a particle in a one dimensional 
random potential. 

We evaluate the path integral ( |L.3|) numerically by mapping it to the associated 
Schrodinger equation. In dimension d = 1 this mapping (see Ref. fl4| Eqs. (3.12)-(3.18)) is 
given by 

-R(L)=R 
'R(0)=R' 

where 



rR(L)=R, 

Z(R,R';L)= [dR(u)]exp(-PH[R(u)]) = (R\exp(-PLH)\R'), (3.1) 

Jr(o)=r> 



H = * d l +^R 2 + V(R). (3.2) 

We compute the matrix element by expanding it in terms of the energy eigenstates of H 

(R\exp(-(3LH)\R') = ^ exp(-f3LE n )$„(/?)*$„ (R'). (3.3) 

n 

In order to compute the eigenvalues and eigenvectors numerically we solve the Schordinger 
equation on a one dimensional lattice of iV sites [ 15| . The lattice Hamiltonian is then an 
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N x N matrix with matrix elements given by 

Hij = ~2M^ OW + <W + (^ 2 (z - N/2f + V(i)) 8 itj (3.4) 

where the lattice spacing is A = S/N, and where S is the system size. Since we are 
interested in the continuum limit A will be kept small. Note that the index i corresponds 
to the position R4 = Ai. The eigenvalues and eigenvectors can now be found directly by 
diagonalizing the matrix using a standard numerical routine ||15|| . Once these are known we 
can construct the partition sum at any value of L using Eq. (ft.3|) . 

The random potential V(R) is generated by first generating a Gaussian correlated random 
potential V^(R) that satisfies 

(V^(R)V^R')) oc exp (-(R - R'f/i 2 ) . (3.5) 

Since we are making a lattice approximation we need a sequence of N numbers {Y^{i)}%=x,..,N 
that obey (V^(i)V^(i + 1)) oc G(l), where in this case G(l) = exp (— A 2 / 2 /£ 2 ). These numbers 
will then be placed on the N lattice sites in the given order. To generate such numbers we 
use a method described in reference [jl~6|] . The procedure is to first generate a sequence of 
N uncorrelated random numbers {U(i)} with a Gaussian distribution. These numbers are 
then fast Fourier transformed, using a standard numerical routine |15| , to yield the sequence 
{U(i)}. Next, we calculate the N numbers defined by W(i) = \J G{i)U{%), where G(i) is 
defined as the Fourier transform of the correlation function G(i). Finally, taking the inverse 
Fourier transform of the sequence {VK(z)}, yields {W / («)}, the sequence with the desired 
correlation function G(i). Now, in order to generate quadratic correlations we choose ^ 
such that the Gaussian correlation function is well approximated by its leading quadratic 
term over the range of the system size. The approximate condition for this to hold is that 
Z/S > 1/V2. In this way we generate a well defined set of random numbers which obey 
approximately the correlation function given in Eq. ( |1.10|) . 

In Fig. |l| we plot a correlation function that is generated by the above method. On the 
same graph we plot the corresponding quadratic approximation. Notice that in this case 
when \R — S/2\ ~ 10 the quadratic approximation begins to deviate from the generated 
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FIG. 1. Plot of (V(S/2)V(R)) versus R. The circles are generated by averaging over 10,000 
samples. The solid line is a plot of the quadratic approximation to the correlation function given 
by Eq. (U). The parameters are £ = 10^6, S = 40, N=200. 

correlation function. This discrepancy turns out to be unimportant as long the quantities 
that are numerically computed (such as end to end distance) do not exceed this range of 
validity. In order to reduce errors due to the finite size of the lattice we found it useful to 
take our sample of random numbers from a set which was about five times N. In all cases 
we tested the reliability of the samples by directly computing the correlation function and 
comparing to the analytical expression for the correlation function given by Eq. ( |3 . 5| ) . 

IV. RESULTS AND DISCUSSION 

We discretized the Schrodinger equation on a lattice of size N = 200. Once the eigen- 
values and eigenvectors are known then we can approximate, for instance the mean squared 
displacement, for each random sample. We then average over the samples to get an ap- 
proximation to the quenched average. For simplicity we set M = 1/2 and (3 = 1 for all 
cases. 

In Fig. |2| we graph the mean squared displacement with one endpoint fixed as a function 
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FIG. 2. Plot of (R^{L)) as a function of L. The parameters are n = 0.3, a = 0.0811, A = 0.2, 
£ = 10\/6. The solid line corresponds to the analytical solution given in Eq. fl2.7|) . The triangles are 
generated by averaging over 1000 samples and the circles represent averaging over 10, 000 samples. 

of L. We do this for two different number of samples in order to check convergence towards 
the corresponding analytical solution. Note that in the labels of the plots the average over 
the disorder is denoted by a second set of brackets rather than an overbar. In Fig. |3| we 



graph (Rt(L)) 2 as a function of L. We use the same parameters as in Fig. |[ It is clear from 
the graphs that the numerical results are consistent with the exact curve. As the number of 
the samples is increased the numerical curves get closer to the analytical solution. In Fig. £| 



we plot [R 2 F (L)) vs. L. This quantity is computed numerically using the expression given 
in Eq. ( |1.5| ). We found that the numerical results were extremely close to the analytical 
prediction after averaging over only 200 samples. 



We now turn our attention to the quantity (Rq(L)), which was discussed in the intro- 



duction. In Fig. j|. we graph (Rq(L)) vs. L. In order to visualize the predicted shift in 
Eq. ( |2.17D we include the exact solution of the zero disorder case. We use the same param- 
eters as the previous figures. Again, we find close agreement between the computational 
results and the analytical solution. The shift due to the disorder is clearly evident and is 
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FIG. 3. Plot of (Rt(L)) as a function of L. The solid line is a graph of the analytical solution 
in Eq. (|2.§|). The triangles are generated by averaging over 1000 samples and the circles represent 
averaging over 10, 000 samples. 




FIG. 4. Plot of (Rp(L)) as a function of L. The solid line corresponds to the analytical solution 
given in Eq. ( |2.10| ). The circles are generated by averaging over 200 samples. All the parameters 
are the same as in Fig 0. 
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FIG. 5. Plot of (Rq(L)) vs. L. The solid line is the analytical solution given in Eq. (ETlTl) . The 
circles are generated by averaging over 8000 samples. The dashed line is the analytical solution for 
zero disorder (a = 0). 



very close to the predicted value. In Fig. | we plot (Rq(L)) 2 vs. L and compare with a 
plot of the analytical solution in Eq. (|2.18 ). For small L there appears to be a discrepancy 
between the data and the analytical solution, whereas for larger L the two curves are very 
close. This is due to the fact that the random potential is generated on a grid with grid 
size of 0.2. Thus for L shorter than 0.2 the particle can not see the random potential and 



(Rq(L)) vs. L should average to zero. Indeed the significant deviation occurs on this length 
scale. 

Finally, we turn our attention to the quantity C(l). It was evaluated numerically from 
the expression given in Eq. ( |1 . 1 1| ) . In Fig. |7| we plot C(l) vs. I with L large and fixed. 

It is clear from the graphs that the numerical solutions are consistent with the exact 
analytical solution. As expected, as the number of samples is increased the numerical results 
get closer and closer to the exact curve. Based on these results we can safely conclude that 
the replica calculation is indeed correct and does describe the averaged properties of the 
polymer. 
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FIG. 6. Plot of (Rq(L)) 2 vs. L. The solid line is the analytical solution given in Eq. ( [2.18 ) 
The circles are generated by averaging over 8000 samples. 
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FIG. 7. Plot of C(Z) vs. I, with L = 10. The solid line is the analytical solution given by 
Eq. ( |2.20| ). The circles are generated by averaging over 250 samples. All parameters are the same 
as the previous graphs. 
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FIG. 8. Plot of (Rt(L)) vs. L. The solid line is the analytical solution given in Eq. (|4.2|). The 
circles are generated by averaging over 1000 samples. We take \x = 0.001, and all other parameters 
are the same as in Fig. ||. 

It is interesting to study the infinite volume limit \i — > 0. Here, the polymer does not 
see the confining harmonic potential and its properties are determined only by the random 
potential. Taking the \i — *■ limit the exact expressions simplify to 



(R£(L)> 



d da 4 



da 



(4.1) 



(4.2) 



Notice that in the no disorder case (a = 0) the later quantity is zero, but once the disorder 
is turned on it scales like L A with a coefficient that is independent of temperature. Also, 
Eq. (|4.1| ) indicates that for small L the polymer wanders diffusively but for larger L it 



wanders much faster than diffusion. In Fig. § we plot (R~t(L)) vs. L when \i is chosen to be 



very small. On the same graph we plot Eq. ( |4.2|) . In Fig. ||we plot (Rj,(L)) — (Rt(L)) 2 vs. 
L. On the same graph we include the analytical prediction. It is clear that the numerical 
results agree well with the analytical predictions. Also, from Eq. (|2 . 1 7]) we see that the 



quantity (Rq(L)) diverges as /i — >• 0. This implies that the boundary conditions on the 
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FIG. 9. Plot of (R^(L)) — (Rt(L)) 2 vs. L. The solid line is the analytical solution given by 
subtracting Eq. ( f4.2| ) from Eq. flhip . The circles generated by averaging over 1000 samples. We 
take fi = 0.001, and all other parameters are the same as in Fig. |2[ 

chain are crucial in determining which quantities are well defined in the infinite volume 
limit. 

The physical consequences of our results are surprising and may seem counter intuitive 
at first glance. The very long range correlations of the random potential lead to a very fast 
wandering of the free end of a tethered chain. However, the deviation from the average 
position (R^(L)) — (R r (L)) 2 does not depend on the random potential. Also, if both ends 
are free to move then the end to end distance (R^(L)) behaves as if there is no random 
potential. This behavior can make more sense if we study the nature of the random potential 
samples that satisfy the quadratic correlation. We find that the typical random potential 
(see Fig. [II]) is smooth and slowly varying on short scales but contains peaks and valleys 
on scales close to the system size. So as L is increased the polymer has a greater tendency 
to be found in the deepest potential well in the sample, which, for large sample sizes, is on 
average located very far from the fixed end of the polymer. So for the case when one end 
is fixed we expect the end to end distance to grow very fast with L, as the bulk (center 
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FIG. 10. The lower curve represents a typical random potential with long range quadratic 
correlations (£ = 10\/6)- For comparison the upper curve represents a random potential with a 
shorter correlation range = 1). 

of mass) of the polymer moves far away from the fixed end. This behavior is the same as 
in the case of short range correlations ||, where the polymer will typically form a tadpole 
conformation with the tail tethered to the origin and the head far away in some region of 
low potential. On the other hand, when both ends are free the entire polymer will simply 
curl up in the region of low potential and the end to end distance should depend only on 
the local behavior of the random potential. Now, since the potential samples are smooth 
and slowly varying on short scales we do not expect the disorder to have much of an effect 
on the local behavior of the polymer. What is very interesting though is the fact that the 
chain that is free to move behaves as if the random potential has no effect at all. Similarly, 
the fluctuations around the average position in the case when one end is fixed turn out to be 
totally independent of disorder. This is expected to be a special feature of the quadratically 
correlated random potential, and is not likely to hold in the general case of long range 
correlations. See discussion below about the crossover to shorter ranged correlations of the 
disorder. 
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It is useful to compare our results with those for directed polymers. Here, the arc-length 
u corresponds to the distance along the directed axis which is a fixed direction in real space, 
and the vector R(w) is the position of the directed line in the transverse hyper-plane. In 
that case the random potential is usually taken to depend on u and satisfy 

(V(R, u) V(R', u')) = S(u - u')f ((R - R) 2 ) , (4.3) 

i. e. the random potential is taken to be uncorrelated along the directed axis, unlike the 
situation described in Eqs. ( |1.1|J1.2| ) were the random potential is independent of u. Parisi 
has shown pi that if / is quadratic, then the mean squared displacement of one end of the 



directed polymer satisfies (Rf,(L)) oc L 3 . This is to be compared with the L 4 dependence 
that we have found for the case when the random potential is independent of u. 

To explain the different scaling we employ a Flory-type argument similar to the argument 



used in |TlJ for directed polymers. Allowing for a rescaling of the arc-length variable u by a 



scale £, i.e. u — > £u and the position variables R(u) — > £^H(£u) we see that the random 



potential which satisfies Eq. ( |1 . 2|) with a correlation function behaving in general like 

/ ((R - R') 2 ) ~ const. ^— (R - R) 2(1 ~ a) , (4.4) 

1 — a 

scales like £ x with 

2A = C2(l-a) (4.5) 

The difference with directed polymers is that in that case one has to subtract a 1 from 
the right hand side of Eq. (|4.5|) because of the delta function in Eq ( |4.3| ). Now in a Flory 
argument one assumes that the two terms in the hamiltonian given in Eq. ( |1.1| ) scale the 
same way (here we consider only the case of \x = 0, since \x ^ breaks scale invariance). 
Since the 'kinetic' energy term scales like £ 2< *~ 2 , and this should be equal to £ x , we see that 

C = ~r~7 • (4.6) 
1 + a 



Thus (R|,(L)) ~ £ 4 /( 1+Q as opposed to L 3 ^ 1+a ^ for directed polymers. In the quadratic 
case a = 0, and we get the L 4 behavior we were looking for. Notice that we derived here 
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a prediction for the behavior of (R|,(L)) for the case of long ranged correlations of the 
disorder which are not quadratic but are characterized by a power law determined by the 
value of a. Thus the power of L will decrease for shorter ranged correlations than quadratic. 
It is interesting to assess the accuracy of the Flory argument in practice. In the function 



(R^(L)) — (R T (L)) 2 or in (R|,(L)) the leading power of L cancells out and one is left with 
a subleading L 1 behavior (for quadratic correlations). 

Within this toy model it is interesting to compare the differences between the annealed 
and the quenched averages. The annealed average applies when the obstacles in the medium 
are randomly placed and mobile. In this case the replica trick is not necessary and the 
random potential can be averaged directly. Alternatively one can use the results obtained 
with Z n but instead of taking n to 0, substituting n — 1. We easily find that 



/ dRR 2 Z(0, R, L) 1 



(R 2 (L)) = J - v I I ' ~ -== , (4.7) 

when L is large, and where we have taken the fi —>■ limit. So in an annealed medium 
with long range quadratic correlations a very long polymer chain will collapse around the 



tethered end. Similarly we find that (R, 2 F (L)) ~ 1/yaL, and so if both ends are free then 
the polymer will collapse in the same way. This behavior is in stark contrast to the quenched 
case where the effects of the random medium is quite different. 

For the case of short ranged correlations and a chain that is free to move one usually 
argues that the annealed and quenched averages coincide in the infinite volume limit. 
This is due to the fact that the system can be divided into subregions, much larger than 
the chain, each containing a different realization of the potential. The moving chain can 
sample all of these and find a realization very similar to the one it induces around itself in 
the annealed case. But this argument does not apply to the case of long ranged correlations 
of the random potential, with the correlation length larger than the system size, since such 
a division to subregions will not yield independent realizations. 

It will be interesting to investigate how the physical properties that we have found change 
as we move towards the regime of short range correlations. In our model we can control 
the correlation length by varying the parameter £ in Eq. ( |3.5|) i.e. in the Gaussian form. 
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For small £ the correlation is certainly not quadratic, and it approaches a 5-function in 
the limit £ — > 0. For arbitrary £ we expect that the average mean displacement squared 



in d — 1 will scale as (Rt(L)) 2 oc L 7 ^. Numerically, we can estimate the exponent 7(£) 



by performing a linear fit to the plot of log (Rt(L)) 2 vs. logL and measuring the slope. 
For short range correlations we found that the numerical method described in Sec. [TTT] was 
unreliable. The reason for this is that the sum over energy eigenf unctions in Eq. (|3.3|) is 
unstable since a typical overlap $ n (i2)*$ n (i?') (for short range correlations) is a number on 
the order of 10~ 15 . However, we were able evaluate Eq. (|3.3|) accurately for all £ by solving 
the Schrodinger equation on a lattice using a fourth order Runge-Kutta algorithm with a 



very small time step (t ~ 10~ 4 ). We found that for large L the quantity (Rt(L)) 2 saturates 
at a constant value due to the finite size of the system, and so we do not expect a power 
law scaling for large L. However, for L sufficiently small (before the onset of saturation) 
the mean displacement squared does obey a power law and a linear fit on a log-log plot was 
excellent for all £. 

In Fig. |11] we plot 7 vs. l/£ 2 for a range of £. For each point we averaged over 8000 
samples on a lattice of size N = 300, and in all cases the strength of the random potential is 
taken to be large (g 3> 1)- We can see from the plot that 7 falls from 4, in the case of very 
long range correlations, to about 2.5 for very short range correlations. The case of delta 
correlated random potentials has been studied by Nattermann [|J using Flory arguments. 
Nattermann finds that for strong disorder (g ^> 1) the mean squared displacement behaves 



like (Rt(L) 2 } oc L 2 ^/g (ln(L)) 3 ^ 2 in d — 1. It is clear from Nattermann's arguments that 



(Rt(L) 2 } ~ (Rt(L)) , and so it is safe to compare our numerical results with his analytical 
expression. So while we find a scaling that is slightly faster than balistic (~ L 2 ' 5 ), Natter- 
mann finds a weakly subballistic behavior (~ L 2 ln(L) ( ~ 3//2 ' ) ). Nevertheless, it is comforting 
to see that both results are fairly close to a balistic scaling (~ L 2 ). 

We now turn our attention to the chain that is free to move. Here, we find that for short 
range correlations (£ < \/5) the end-to-end distance rises linearly for small L and saturates 
at a constant value for large L. This saturation is not due to the finite size of the system 
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FIG. 11. Plot of 7 vs. l/£ 2 . 

since it occurs at a value of L far less than the length at which a free chain would saturate. 
Typically, for L > 1 we find that (R F (L)) oc L° as compared to (R F (L)) oc L when the 
correlations are long range and quadratic. In order to quantify this crossover between the 
long and short range behavior we assume that the scaling relation (R F (L)) oc L 5 ^ holds 
for L > 1. Again we can estimate the exponent by measuring the slope of the line in a 
linear fit of log (R F (L)) vs. logL. 

In Fig. [12] we plot 5 vs. l/£ 2 for a range of £. We computed the end-to-end distance for 
L in the range 5 < L < 10. For each point we averaged over 1000 samples on a lattice of size 
N = 200, and in all cases the strength of the random potential is taken to be large (g 3> 1). 
We can see from the graph that as the correlation range is decreased S falls rapidly from 
about 1 to a value close to zero. This implies that the behavior of (R F (L)) is very strongly 
dependent on the correlation range. These results are consistent with the Flory arguments 
in Ref. |3|,|]] where it is predicted that a long polymer in a delta correlated random potential 
will have fixed size i.e. R? ~ L° in a sample of finite volume, and with the variational results 
of Ref. § (see also 
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FIG. 12. Plot of 5 vs. l/i 2 



V. CONCLUDING REMARKS 



In this paper we presented a model of a polymer chain in a quenched random media which 
was exactly solvable using the replica method. The analytical results were subsequently 
found to be in close agreement with a numerical solution in d = 1 . Based on these results we 
can safely conclude that the replica method is accurate in describing the averaged properties 
of the polymer. The physical picture that emerged was interesting and somewhat surprising. 
We found that a quadratically correlated disorder has a major effect on the size of a polymer 
with one end fixed, but has no effect on the size of a chain that is free to move and find an 
optimal position. We also found that the quenched and annealed cases are rather different: 
In the annealed case a long chain collapses to a point. 

Overall, we have learned that chain properties depend strongly on the correlation range 
of the random media. However, there are still some open problems. For instance, it would 



be very useful to have an analytical derivation of the results depicted in Figs. |TT| , [12| . Also, it 
may be fruitful to investigate various non-equilibrium properties of polymer chains in long 
range correlated random media, such as transport properties and chain dynamics. We hope 
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that our results for the simple case of quadratic correlations will be a useful starting point 
for a more detailed analysis of these problems. 
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APPENDIX: 



Here we show some of the intermediate steps that lead to Eq. (|2.7| ). We first write the 
replicated partition sum as 

Z n ({R a }, {R'J; L) = J d\e- X2 / 2 Z{{K a }, {R'J; L, A), (Al) 



where 



Z({Rj,{Rj;L,A) = H / [dRa]e- 

a=1 JRa(0)=R a 



l3H a (X) 



After performing the path integrals using equations (3.39 — 3.41) in [14], we get 



Z({R a },{R' a };L,X)=N e-^, 



(A2) 



(A3) 



where 



with 



A = ^^/ZMcoth (Ly/fi'/M 



B = ^/JTM ^cosh \Ly/fi'/M) - lj ^sinh [Ly/n'/M 
1 



D = B - L///2, 



(sinh (L^n'/M 



a = — -Jo A. 



(A4) 



(A5) 
(A6) 
(A7) 
(A8) 
(A9) 



The exact form of the normalization Nq is unimportant as it will cancel out later. The next 
step is to perform the Gaussian integrals over the A variables in Eq. (M). This yields 



Z n ({R a }, {Rj; L) = Ni exp |-C/ (E R ' + E R ?) " V (E R ° + E ^ 

-2H/^R a -Rl}, (A10) 



with 



25 



U = /3A 



V 



H' 2 + Sf3anD 

w = pc. 

Now, for arbitrary n we can write 

/ dRf-dR n R 2 Z n {{0},{R a };L) 



(All) 
(A12) 
(A13) 



(Rf(L)> 



/ dR 1 ---dR n Z n ({0},{R a };L) 
~ndU ln I dR i--- dR nZ n ({0},{R a } ] L) 



Id d 
~n2dU 11 



7T" 



d ( 1 
2 



,[/ U(U + Vn)J' (AM) 
We have used the fact that the eigenvalues of the matrix associated with the quadratic 

form in Eq. ( |A10|) are U with multiplicity n — 1 and Z7 + nV with multiplicity 1, and the 

determinant is the product of the eigenvalues. The n — > can now be safely taken to yield 

(A15) 



(R^(L)> = lim(R;> = 



d ad f B 
+ — ■ 



.4 



which simplifies to yield Eq. 



To calculate (Rr(L)) we use 

^ JdR 1 ---dR n R 1 -R 2 Z n ({0},{R a };L) 



(R T (L)Y 



J dR 1 ---dR n Z n ({0},{R a };L) 



d d 
n(n- 1) ~~ dV 
1 d ( d d 



In / dRi • • • dRnZ n ({0}, {R a }; L) 



n(n - 1) 2 \dU dV 



In 



7T 



U n - l (U + nV) 



V 



U(U + Vn) 



(A16) 



Next, we show how to calculate (Rf,(L)). 

/ ]1 dR a n ^ (Ri - R'i) 2 Z n ({R a }, {R' a }; L) 



f n dR a n dR'„^n({Ra}, {R'J; 

In y J] rfR a J] rfRlZ n ({R a }, {R^}; L) 



n \dU dW 
Id f d d 



n2 \dU dW 
d 



In 



7T 



2n 



(U + W)"" 1 ^ - W) n (U + W + 2nV) 



d 



U-W P(A-C) 



(A17) 
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which yields Eq. (|2~T0|) . 



Finally, we calculate (Rq(L)). This is given by 

J dUi ■ ■ ■ <iR„R 2 Z n ([R a }, {R a }', L) 



J dRi ■ ■ ■ dR n Z n ({R a }, {R a }; L) 
— -jjj In / dR 1 ---dR n Z n {{R a },{R a };L) 



1 d d 
~2^2dU 
d ( 1 



In 



7T 



[2(U + W^-^U + W + 2nV)_ 
2V 



4\U + W (U + W){U + W + 2nV) 
In the limit n — > we obtain 



d 



2V 



A\U + W (U + W) 2 

2 



d a 



B 



/i 2 \A + C 



d 



A(3(A + Cy 



(A18) 



(A19) 



which gives rise to Eq. ( p.l7|) . (Rq(L)) 2 is calculated similarly from (R x • R 2 ) and one 
obtains 4a d/[i 2 in agreement with Eq. (|2.18 ). 
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